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£h ' Abstract 
£3 ! 

The structure of 8 He is investigated within a three-cluster microscopic model. The 
three-cluster configuration a + 2 n + 2 n was used to describe the properties of the ground 
state of the nucleus. The obtained results evidently indicate the existence of a neutron 
halo in 8 He. 

t-H ! 

O 

1 Introduction 

o. 

The development of the experimental technique made it possible to investigate light nuclei 
with large neutron excess, i.e., the nuclei for which the ratio 7] = (N — Z) /A is significantly 
larger than for common ones. Such nuclei lie near the drip line and are /3-unstable. They 
live a short time and transform by emitting electrons into nuclei with approximately equal 
number of protons and neutrons. A number of unexpected properties were discovered 
in those nuclei, for instance, a neutron halo. It is natural that many attempts were 
undertook to explain those properties within microscopic and semi-microscopic methods 
0,gl,|,|,|J,l. 

Our aim is to investigate the structure of the 8 He ground state. It is interesting that 
8 He has the largest value of r] = 0.5 among other nucleon-stable nuclei. Note that the 
average values of r/ = 0.4 for nuclei near the neutron drip line. As early as in 1960, 
Ya. B. Zeldovich || and V. I. Goldansky indicated a possibility of the existence of 
8 He isotope. It was experimentally confirmed [11| in the middle of sixties. The subsequent 
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analysis shows that the lowest threshold of He decay is He+2n and lies 2.1 MeV above 
the ground state and energy of the threshold a + 4n equals 3.10 MeV (see, for instance, 

The most complete information on light nuclei with neutron excess can obtained with a 
microscopic model. In this case, the problem is connected with solving the many-particle 
Schrodinger equation with a fixed (chosen) nucleon-nucleon interaction. The equation has 
to be solved with some simplification based on one or other physical considerations. The 
Resonating Group Method |l3|] or its Algebraic Version [14, |l5| is one of such methods. 



In this paper, we make use of the Algebraic Version of the Resonating Group Method 
in which 8 He is considered as a three-cluster configuration a + 2 n + 2 n. It is obvious that 
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we make a priori some assumptions on the structure of the nucleus. First, wave functions 
of each cluster are modelled by the shell-model functions. Second, valent neutrons unite 
in dineutron clusters. As a justification for such an assumption can be served the fact, 
indicated by A. B. Migdal |jiq 1, that the interaction between two neutrons may be increased 
significantly in the presence of the third particle. It can give rise to the creation of 
the dineutron clusters on the surface of a nucleus. The chosen clusterization allow us 
to consider an a-particle as a core, despite that the lowest threshold of 8 He decay is 
6 He+n+n. Earlier, A. A. Ogloblin IlTj indicated the importance of a cluster configuration 
a + An. He pointed out that the bound energy of two neutrons in 8 He is two times larger 
than that in 6 He. This fact led him to the conclusion that 6 He cannot serve as a core and 
the neutron halo in 8 He has to be consisted of four neutrons. 

Note also that the usage of dineutron clusters is a quite grounded approximation. For 
example, in fll^l, dineutron and also diproton clusters were successfully used to describe 
exit channels^ the reactions »H+»H J**H-n + n and ^Re^Re+p+p respectively. 
Besides, in [19], main features of n Li was reproduced within the cluster configuration 
9 Li+ 2 n with a pointless dineutron. 



2 Method 

The present method for investigation of the 8 He ground state is based on the Algebraic 
Version of the Resonating Group Method (AV RGM). For a long time, this version was 
used for studying the bound states of two-cluster systems, reactions with a few open 
channels, interaction of these channels with collective monopole and quadrupole modes, 



and also processes of full disintegration of light nuclei [20, 21, 22]. Recently, the AV RGM 
was actively applied to describe three-cluster systems [^, 24, 25, |26[| . 

The Algebraic Version of the Resonating Group Method is based on the usage of an 
oscillator basis for solving bound state problems and problems of continuous spectrum 
states. This is achieved by expanding a wave function of inter-cluster motion in the 
oscillator basis. As a result, a trial three-cluster function takes the form 

*(A) = £ C n A [$ 1 (A 1 )$ 2 (A 2 )$3(A 3 )/n(qi ! q2)] , (1) 

n 

where A is the antisymmetrization operator, <£j(Aj) are the internal functions of the 
cluster, which are selected in one or other form prior to solving the problem (for instance, 
in the form of many-particle oscillator shell functions as in our case) ; the set of coefficients 
C n is nothing else but a wave function in the oscillator representation. This function 
should be obtained from a system of linear equations: 



< n, | H | ri > -E < n \ ri > C n > = 0, (2) 

which is derived directly from the many-particle Schrodinger equation. The oscillator 
functions /«(qi,q2), where qi and q 2 are Jacobi vectors fixing a position of clusters in 
space, are determinate in the six-dimensional space and constitute the irreducible repre- 
sentation [iVOO] of the unitary group U (6). Thus, the composite index n consists of indices 
(six in total) of the irreducible representation of the U(6) group and its subgroups. 
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The choice of one or other reductions of the U(6) group is dictated by considerations of 
physical lucidity and simplicity of numerical realizations as well. To consider the bound 
state problem, it is convenient to use bases, whose classification is connected with the 
following reduction of the U(6) group: 

U(G) D U(3) ® U(3) \Nih,N 2 ,h,LM> 

u u 

50(3) ® 50(3) D 50(3) 
17(6) D 5C/(3) (8) 17(2) => | (Xfx)u, ujLM > 

u u 

50(3) 0(3) 

The first basis is usually called the basis of two uncoupled oscillators or bioscillator 
basis (BO). Each of the SU(3) groups, associated with one of the Jacobi vectors qi and 
q2, generates the quantum numbers N\, l\ and N 2 , l 2 . They are the principal quantum 
number (or the number of oscillator quanta) and partial angular momentum along the 
respective Jacobi vector: 

\Ni,h,N 2 ,l 2 ;LM> 

The second basis is an "577(3)" basis. Wave functions of this basis are classified 
through the well-known Elliott indices (A,//) of the SU(3) group, multiplicity index u 
arising in the reduction SU (3) C 50(3), and quantum number v = 2 (Ni — N 2 ) connected 
with oscillator quanta along the Jacobi vectors (fiand q 2 : 

\(Xfj,)f,uLM > 

The total number of oscillator quanta equals N = N\ + N 2 = A + 2/i and defines the 
irreducible representation of the U(6) group. For a given N, i.e., for a fixed oscillator shell, 
functions of both bases related to each other through a unitary transformation, because 
these bases are eigenfunctions of the same oscillator hamiltonian in the six-dimensional 
space. Thus, they are equivalent. Note that the unitary matrix connecting these two bases 
consists of the Clebsch-Gordan coefficients of the SU(3) group for the decomposition of 
the product (iViO) ® (JV 2 0) (A/i). Thus 

(A^)z^; ujLM >=J2 U (Ni,h,N 2 , l 2 ; (X/jl)u;u) \N 1} h,N 2 , l 2 ;LM > 
h,h 

However we make use of two bases. This is because the bioscillator basis has more 
natural quantum numbers. Meanwhile, the SU(3) basis is more convenient for numerical 
implementation, in particular, for eliminating Pauli-forbidden states. Besides, the usage 
of two bases gives additional information on optimal subspaces, which allow one to obtain 
reliable results with minimal effort. 

The elimination of Pauli-forbidden states is performed by diagonalization of the matrix 
of the antisymmetrization operator 

\\<n\n>\\, (3) 

calculated between the basis functions (||). Pauli-forbidden states correspond to those 
eigenfunctions of the matrix ||<n|n'>[| = ||<n A n' > || which have zero eigenvalues. 
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Pauli-allowed states are a combination of original basis functions of a given oscillator shell 
which are eigenfunctions of the antisymmetrization operator. It should be noted that the 
matrix || < n \ n' > [| has a block structure. Non-zero matrix elements correspond 
to overlapping basis functions of the same oscillator shell, i.e. those oscillator functions 
which obey the condition N = N' . To solve the Schrodinger equation in matrix form, one 
has to eliminate Pauli-forbidden states. Let us e a and {U"} be respectively eigenvalue 
and eigenfunction of the antisymmetrization operator. Then, the system of equations (^) 
should be transformed to the representation of Pauli-allowed states: 



E 

a' 



< a \ H 



o 



> -E5 n 



C a > = 0, 



(4) 



where 



< a I H I a 1 > 



is a matrix of hamiltonian between Pauli-allowed states connected 



with the matrix 1 1 < n H n' > 1 1 by the relation 



H 



U° 



n 



H 



n 



U5 



In this connection, for the bioscillator basis the original scheme of classification is 
totally changed, but the quantum numbers (A//) are preserved for the SU (3) basis, because 
the matrix 1 1 < n \ n' > 1 1 is off-diagonal with respect to the quantum number v only. 

We omit all details of matrix elements calculations of the microscopic hamiltonian and 
antisymmetrization operator, by referring reader to the paper [27] where one can find 
basic formulae and recurrence relations for matrix elements of operators of the physical 
importance between bioscillator functions. 



3 Results 



Results, represented in this chapter, were obtained with the Volkov potential [28|. The 
only free parameter, oscillator radius ro, was chosen to minimize the threshold energy 
of the 8 He decay into He and two dineutrons. It turns out to be 1.51 fm. Under such 
conditions, the energy of the 4 He+ 2 n + 2 n threshold equals -22.15 MeV and the bound 
state energy of ft-particle is -26.84 MeV. Coulomb interaction was neglected because it 
leads to a shift of the bound state energy and threshold energy by the same value. 

In what follows, we use two different trees of Jacobi vectors. In the first tree which 
we call "T"'-tree, the vector qi defines the distance between two dineutrons, and the 
vector q2 fixes the distance between the center of mass of two dineutrons and a-particle. 
The second tree is called Y-tree. In this tree, the first vector qi determines the distance 
between a-particle and one of the dineutrons, and the second vector q2 is connected with 
the distance between the second dineutron and center of mass of the first dineutron and 
ft-particle. 

As we concern with the ground state only, then we need to use wave function of the 
S'-state. In this case oscillator basis is reduced significantly. For instance, the bioscillator 
basis involves oscillator functions with even values of N\ and A^- Besides, partial angular 
momentum l\ = l<i- Actually, we need only three quantum numbers to classify basis 
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Figure 1: 8 He ground state energy as a function of the number N of oscillator shells involved 
in the calculation. 

functions with L = 0. They are (N±, N2,l = h = h) for bioscillator basis, and (A/i, v) for 
SU(3) basis. 

Bound state energy and optimal subspaces. The ground state of 8 He is con- 
sidered with the basis which involves all oscillator functions of 15 lowest oscillator shells, 
i.e., basis functions with even values of the principal quantum number N up to N = 30. 
The total number of the original basis functions equals 815 and the total number of Pauli- 
allowed states reduces to 399 functions. Such a number of basis functions provides a fairly 
good convergence of the bound state energy, as is demonstrated in Fig. |l[ In this figure, 
we display the ground state energy as a function of the principal quantum number N . The 
energy is counted from the threshold a + 2 n + 2 n. In Fig. [I], we also display the ground 
state energy obtained with some subspaces of the total space of the oscillator basis used. 
In the BO basis, such a subspace is defined by the maximal value of the partial angular 
momentum I = 0, while, for SU (3)-basis, such a subspace involves basis functions with 
H < 4. The later subspace consisting of 274 functions gives the energy which noticeable 
differs from "exact" one, obtained with the total basis. But, with the former subspace in- 
cluding only 118 functions, we obtain the energy which is very close to the "exact" value. 
This is probably connected with that the interaction between clusters is most strong in 
the S-state. 

In Fig. ^, we display the wave function of the 8 He ground state, more exactly, the 
coefficients C a of expansion over Pauli-allowed states. Two labels N s h and N a . s are used to 
classify Pauli-allowed functions (a = {N s h, N a . s }). The first label N s h numerates oscillator 
shells and the second one, N a _ s , , numerates Pauli-allowed states of a given oscillator shell. 
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Figure 2: Wave function of the 8 He ground state in oscillator representation. 
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Figure 3: Contribution of different oscillator shells N to the wave function of the 8 He ground 
state. 



The expansion coefficients C a were determined in the 5Z7(3)-basis, where the indices (Xfj,) 
are good quantum numbers after eliminating Pauli- forbidden states. The detailed analysis 
shows that the main contribution (around 80%) to the wave function comes from the basis 
states with [i = 2, while the basis states with /x = give only 9%. Note that the former 



shells are displayed) that the main contribution comes from the lowest oscillator shells, 
however the contribution of shells with large N is also noticeable. It indicates a substantial 
clusterization of the nuclei, i.e., for a large amount of time, valent neutrons move far from 
the a-particle, making a neutron halo. 

In order to obtain additional information on the role of different subspaces of the total 
space of oscillator functions, we impose various restrictions on the quantum numbers of 
basis states. First, for the bioscillator basis we took a subspace with the maximal value of 
partial angular momenta (I = l\ = l%) I = 0, I = 2 and I = 4. It was made for both Y- and 
T-trees of Jacobi vectors. For the SU (3)-basis, we used only T-tree and the restriction 
was imposed on the maximal value of \x = 0, 2 and 4. Results of such calculations are 
presented in Table [j]. One can see that subspace l\ = I2 < 2 for y-tree of the bioscillator 
basis is the most optimal part of the total basis, because 54% of the total basis (or 219 
functions) gives the ground state energy very close to the "exact" value. 

Effects of the Pauli principle. To understand the role of the Pauli principle in a 
three-cluster system, we investigate the contribution of different Pauli-allowed states to the 
wave function of ground state. Each of Pauli-allowed states, being an eigenfunction of the 




than 93%. 



It is seen from Fig FJ and more clearly from Fig. || (where weights of different oscillator 
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Table 1: Ground state energy of 8 He counted from the threshold 4 He + 2 n + 2 n. 



T$ £1*31*3 

LJ OjQ-LQ 


Taroni 


i .J I.I IJCIj^C 


E, MeV 


TVn lnhpr 




tree 






of functions 


BO 


T 


Xotal 


-2.065 


399 






Z = 


— 1.832 


118 






I < 2 


—2.047 


219 






1 < 4 


-2.064 


293 


BO 


Y 


Total 


-2.065 


399 






Z = 


-1.839 


118 






Z < 2 


-2.065 


219 


SU(3) 


T 


Total 


-2.065 


399 






li = 


6.341 


91 






H<2 


0.172 


196 








-1.335 


274 



antisymmetrization operator, can be marked (characterized) by corresponding eigenvalue 
of this operator. As we mentioned above, the antisymmetrization operator overlaps only 
those functions which obey the relation N\ + N 2 = N[ + N' 2 , i.e., basis functions of the 
same oscillator shell. Analysis of the eigenfunctions shows, that the diagonalization of the 
matrix 1 1 < n A n' > 1 1 reveals states with definite eigenvalues of the antisymmetrization 
operator for a two-cluster subsystem. This means, in particular, that the Pauli-forbidden 
state of three-cluster system is a state when at least one pair of clusters is in Pauli- 
forbidden state. For example, for the two-cluster subsystem a+ 2 n, the oscillator functions 
with the number of oscillator quanta (along the inter-cluster coordinate) N = and N = 1 
are Pauli-forbidden states. In the subsystem 2 n + 2 n, where symmetry of the subsystem 
allows only even functions, we have only one forbidden state with N = 0. As for Pauli- 
allowed states for three-cluster system, they describe the states of the system, when all 
pairs of two-cluster subsystems are out of the Pauli-forbidden region. 

To prove these statements, we consider eigenvalues of the antisymmetrization operator 
(which we denote by \ a ) for Pauli- allowed states and expansion coefficients C a over these 
states for oscillator shell N s h = 20. These quantities, obtained in the SU (3) basis, are 
displayed in Fig. |I[ Seven first functions correspond to the SU (3) irreducible represen- 
tation (A/i) = (20,0), next eight functions belong to the SU(3) irreducible representation 
(A/u) = (16,2) and so on. Last function has the SU(3) symmetry (A/u) = (0, 10). In this 
Figure by dashed horizontal lines we indicate the eigenvalues of the antisymmetrization 
operator for two-cluster subsystem a + 2 n. Note that corresponding values for the sub- 
system 2 n + 2 n equal 1. One can see indeed, that some eigenvalues of the operator A 
for three-cluster system coincide with the eigenvalues of this operator for the subsystem 
a + 2 n. Besides, one notices that such states, corresponding to the even oscillator quanta 
N = 2 and N = 4 in the a + 2 n subsystem, play dominant role in the ground state of 8 He. 

RMS radii. In Table we compare the calculated mass, neutron and proton root- 
mean-square (RMS) radii with available experimental data. The theoretical values of 



Figure 4: Eigenvalues of the antisymmetrization operator Ajv . s . and expansion coefficients C a 
for the oscillator shell N = 20. 
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Table 2: Root mean square radii for the ground state of s He. 



RMS, fm 


Theory 


Experiment ||29J 


Experiment ||30|| 


RMSm 


2.73 


2.37 ±0.18 


2.52 ±0.03 


RMSp 


2.08 


1.89 ±0.17 


2.15 ±0.02 


RMSn 


2.91 


2.50 ±0.19 


2.64 ±0.03 


RMSn-RMSp 


0.84 


0.61 


0.49 



Table 3: Root mean square radii of He, obtained by different methods. 



RMS, fm 


AV RGM 


RRGM [1| 


Shell Model § 


RMSm 
RMSp 
RMSn 


2.73 
2.08 
2.91 


2.41 
1.71 


1.684 



RMS radii are a little larger than experimental ones which we took from pi and pfj. 



This is perhaps because the calculated binding energy is a little less than the experimental 
one. But the present model correctly reproduces the general picture of 8 He. One sees 
that radius of neutron matter is larger then the one of proton matter. The difference of 
these radii is 0.84 fm. These results also indicate the existence of a neutron halo in the 
nucleus. 

In Table g, we collected the mass, neutron and proton RMS radii, obtained with 
different theoretical methods: AV RGM (present calculations), the Refined Resonating 
Group Method (RRGM)|J, and Multi-configuration Shell Model f§. 

Shape of three-cluster system. Having calculated the coefficients C n , a wave 
function in the oscillator representation, we thus obtain the wave function for relative 
motion of the three-cluster system: 

$(qi,<l2) =J2C n (f) n (q, 1 ,ci2). (5) 

n 

By using (||), we can evaluate mean distance between clusters. For this aim, one has 
to calculate the following quantities: 



Ql = J dqidq 2 $*(qi,q2) q? $(qi,q 2 ) 
Ql = J d qi dq 2 $*(qi,q 2 ) q 2 $(qi,q 2 ) 



which, within the regard for the normalization of the wave function and definition of 
Jacobi coordinates, define the sought parameters. For instance, for T-tree, the mean 
value of q\ is connected with the base of an isosceles triangle and the mean value of g 2 is 
connected with its height. 
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Figure 5: Shape of the triangles which are composed by three clusters in He and He. 

The mean distance between two dineutrons turns out to be 2.33 fm and the mean 
distance between a-particle and the center of mass of two dineutrons is 1.42 fm. Thus, 
in 8 He, three clusters form an isosceles, almost rectangular triangle with a-particle at the 
vertex of the right angle. 

Note, that the situation is somewhat different for 6 He. Clusters form an acute-angled 
triangle. Two valent neutrons in the presence of a-particle make a subsystem with the 
RMS radius equal to 2.52 fm which is less than the RMS radius of a free deuteron (2.69 
fm) calculated with the same potential and the same number of basis functions. These 
triangles are displayed in Fig. || 

The difference in the geometry of cluster's disposition in 6 He and 8 He is more likely 
connected with the Pauli principle. There is an effective repulsion between two dineutrons, 
arising from the Pauli principle, which strives to place dineutrons on different sides the 
a-particle. Contrary to the case of 8 He, valent neutrons with opposite orientations of 
spins may unite in a rather compact subsystem in 6 He (due to the presence of a-particle) . 

Density distribution. Proton, neutron and mass density distributions also confirm 
the existence of a neutron halo in 8 He. As seen from Fig. ^, where we display the proton, 
neutron and mass 

density distributions for both 8 He and 6 He, the size of a neutron cloud is substantially 
larger than the size of a proton cloud. Besides, main part of neutrons in 8 He move on the 
surface of the nucleus. One sees the depression in the neutron density distribution at small 
values of the coordinate r. This is due to the Pauli principle, which makes four neutrons 
(united in two dineutrons) move at a relatively large distance from the a-particle. 
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Figure 6: Proton, neutron and mass density distributions in 8 He and 6 He. 
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4 Conclusion 



In this paper, we have investigated the ground state properties of He within the three- 
cluster microscopic model. The three-cluster configuration 4 He+ 2 ra + 2 n was used to 
simulate the dynamics of the eight-nucleon system. The model suggested describes rea- 
sonably well parameters of the ground state: binding energy, mass, proton and neutron 
root-mean-square radii. The analysis of the system shows, that valent neutrons move at 
a large distance from a-particle, forming a neutron halo in 8 He. 
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